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ABSTRACT 

The structure of accretion curtain flows, which form due to disc disruption by a 
strongly magnetic star, is considered. It is shown that a sub-Alfvenic, magnetically 
channelled flow is consistent with matching the magnetic field across the curtain base 
where it meets the disrupted inner region of the disc. The resulting angular velocity 
distribution in the curtain flow is calculated, together with the consequent angular 
momentum tranfer rate to the star. 

It is shown that the transition of material from the diffusive disc flow to the 
channelled curtain flow results in some angular momentum being fed back into the 
disc. This is consistent with the total angular momentum balance, and can result 
in a significantly smaller accretion torque acting on the star than that given by the 
standard model which assumes no angular momentum feed-back to the disc. The sonic 
point coordinates are found and a critical stellar rotation rate results below which the 
sonic point merges with the curtain flow base, due to a reduced centrifugal barrier to 
the flow. Hence, at these lower rotation rates, little thermal assistance is required for 
matter to make the transition to a supersonic accretion flow. 
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1 INTRODUCTION 

1.1 Types of systems 

An accretion curtain forms when a strongly magnetic star disrupts the inner part of its surrounding accretion disc. Material 
leaves the disc and becomes channelled by the stellar magnetic field until it reaches the star's surface, via localised columns 
of shock-compressed gas. The angular momentum transferred in the curtain flow plays an important role in determining the 
spin evolution of the star. Such accretion is believed to occur onto a protostar in T Tauri systems (e.g. Bouvier et al 2007), 
a neutron star in X-ray binary pulsars (e.g. Campbell 1997) and a white dwarf in the intermediate polars (e.g. Warner 1995; 
Frank, King & Raine 2002). 

The intensity of radiations emitted by the shocked hot gas in the accretion columns is rotationally modulated, and 
observations of this yield information on the spin behaviour of the star (e.g. Chakrabarty 1997). Spin- up and spin-down are 
observed, as well as more erratic period changes (e.g. Bildsten et al 1997). The nature of the curtain flow and the associated 
angular momentum transfer are of fundamental significance in these important systems. 

1.2 Previous Work 

Pringle & Rees (1972) suggested that a strongly magnetic neutron star accreting via a disc is likely to disrupt the disc far 
from the stellar surface inside a corotation radius, zuco, at which the Keplerian angular velocity of material in the central 
plane of the disc equals the stellar angular velocity. They argued that inside the disruption radius ionized matter would be 
channelled by the strong magnetic field of the star before being accreted on to its surface. 

Ghosh & Lamb (1978, 1979a, b) constructed a basic model to describe the interaction of a magnetic neutron star with its 
accretion disc. They pointed out that the stellar magnetic field is likely to penetrate the thin turbulent disc as a consequence of 
instabilities and reconnection. The magnetosphere is expected to nearly corotate with the star, and the vertical shear between 
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this surrounding region and the disc material leads to the generation of toroidal magnetic field firom the poloidal component. 
The resulting magnetic stresses lead to angular momentum exchange between the disc and star. The toroidal magnetic field 
reverses sign across vjco, and hence the direction of angular momentum exchange between the disc and star also reverses. 
Outside Wco the star adds angular momentum to the disc, while inside this radius angular momentum is removed from the 
disc. Hence these regions make spin-down and spin-up contributions, respectively, to the stellar torque. 

It was sugcssted by Ghosh & Lamb that magnetic stresses would result in disc disruption in the region of a spherical 
Alfven radius, vja, inside Wco. It was argued that disruption should occur over a narrow annulus, through which material 
leaves the disc vertically and enters a channelled fiow guided by the relatively strong stellar magnetic field. However, no 
mechanism was given to explain disruption of the disc or the vortical transition of the disc flow into the field-controlled 
accretion region. The angular velocity of disc material was taken to change from a Keplerian value to the stellar value across 
the inner boundary layer in which disruption occurs. All the angular momentum of material leaving this layer was assumed 
to be transferred to the star via the channelled flow and a consequent accretion torque was calculated. A stellar spin rate can 
be found at which the magnetic and accretion torques balance, so it is possible to have zero total torque acting on the star 
while it is accreting material from the disrupted disc. 

Campbell (1992) considered the inner part of the disc, inside zuco- It was shown that the vertical equilibrium cannot be 
maintained if the magnetic torque signiflcantly exceeds the viscous torque. This work indicated that disc disruption is due to 
a thermal process which causes the disc to expand vertically over a narrow radial region. 

Most authors have assumed, like Ghosh & Lamb, that all the angular momentum of the magnetically channelled material is 
transferred to the star (e.g. Wang 1987; Armitage & Clarke 1996; Matt & Pudritz 2005; Tessema & Torkelsson 2010; CampbeU 
2011). However, Shu et al (1994a,b) and Ostriker & Shu (1995) developed a modified model of the star-disc interaction in 
which the accretion torque on the star is greatly reduced. In this model the field lines thread the disc, but disruption occurs 
close to the corotation radius. As in the Ghosh & Lamb model, matter is channelled onto the star through a magnetic field 
guided fiow orginating from a thin disruption region in the disc. However, a wind fiow is taken to develop from the disc 
just beyond this region and this carries a significant amount of angular momentum away from the system. Consequently the 
accretion torque delivered to the star is significantly reduced. The accretion and wind fiows emanate from an 'X-type' region 
centred on Wm-, in which it is assumed that a steady state can be maintained. This model was constructed for particular 
application to protostellar systems (T Tauri stars). 

Li, Wickramasinghe & Ruediger (1996) argued that the accretion torque could be effectively zero due to all the angular 
momentum of the field channelled flow being fed back into the disc. They argued that the transfer of angular momentum via 
the channelled flow would require large values of the flold winding ratio \B^/ B^\ at the stellar surface, which would lead to 
instabilities. However, Wang (1997) pointed out that only small values of this ratio are required at the star's surface since the 
torque is related to the stress B^Bp which is sufficiently large even though \B^/Bp\ is small. 

Numerical simulations of the field channelled curtain fiow have been performed, with one of the most extensive by 
Romanova et al (2002). These simulations showed that funnel fiows develop in the inner disc region. Material becomes 
channelled by the stellar magnetic field and this leads to a spin-up torque on the star. The case of a dipolar field was considered, 
with the magnetic moment aligned with the spin axis of the star. Cases with tilted dipoles and multipole components were 
analysed numerically by Romanova et al (2003). These calculations were done in the context of T Tauri stars, in which 
observations suggest that the magnetic field is more complex than dipolar (e.g. Jardine et al 2002). 

Erkut & Alpar (2004) derived a form for the dynamical viscosity from the steady angular momentum balance for a 
Keplerian disc. They then used this in the angular momentum equation to estimate deviations of f2 from inside Wco. For 
sufficiently slow values of the stellar rotation rate, fi*, they found that f2 can change from Qk to f2, over a broader region 
than the boundary layer suggested by Ghosh & Lamb. However, for larger values of fl» boundary layer behaviour is exhibited. 
A detailed disc solution was not found, nor a mechanism for disc disruption. 

It was shown in Campbell (2010) that disc disruption occurs inside the corotation radius as a result of increasing magnetic 
heating. The consequent growing central plane pressure cannot bo balanced by the vertical component of the stellar gravity 
and the disc expands over a narrow radial region. At disruption the magnetic torque is rapidly increasing and material is 
diverted from the disc into a magnetically channelled region, forming a curtain-type ffow. The increased vertical pressure 
gradients allow material to surmount the effective potential barrier and flow through the sonic point. Beyond this point, the 
component of the stellar gravity along the poloidal magnetic held exceeds that due to the centrifugal force and material is 
subsequently accelerated to speeds approaching free-fall values. An approximate expression was found for the position of the 
sonic point, and the energy required to reach it. 



1.3 Paper outline 

The present paper looks at the nature of the accretion fiow in more detail than previous studies. In Section 2 the magnetic and 
dynamical problems are addressed. The equations are solved for the flow and magnetic field, as a function of the stellar rotation 
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rate. Matching of the disc and curtain structures is shown to have fundamental consequences for the angular momentum 
transport to the star. 

In Section 3 the slow magnetosonic point is considered and its coordinates arc calculated. This is shown to be essentially 
the same as the sonic point. A critical stellar rotation rate results at which the slow point approaches the disc surface, 
hence facilitating the transition to supersonic flow. The resulting accretion torque is calculated in Section 4, as a function 
of the stellar rotation rate. This torque can be significantly lower than that given by the standard form, since some angular 
momentum can be transferred back to the disc in the transition layer at its inner edge. Section 5 discusses the results. 



2 THE GOVERNING EQUATIONS 

A dipolar magnetic field is considered since, in the region of interest, this will be the dominant component in white dwarf and 
neutron star systems. In general, the magnetic moment will be inclined to the rotation axis but, as illustrated in Campbell 
(2011), the magnetic field can be split into timo-depondont and time-independent parts with the latter being axisymmetric. 
The time- independent field dominates in the relevant part of the accretion flow and, in cylindrical coordinates (ro, «), its 
unperturbed components are 

where Bo and R are the polar magnetic field strength and the stellar radius, respectively. It will be shown that the accretion 
flow only causes small distortions of the stellar magnetic field from this form. 

The steady magnetohydrodynamic equations will be employed for the curtain flow, in the limit of high magnetic Reynolds 
number. The quantities conserved along poloidal magnetic field lines, that are presented in this section, are similar to those 
occurrring in a steady, axisymmetric magnetic wind flow (e.g. Campbell 1997, Mestel 1999 and references therein). However, 
some important differences arise since the poloidal velocity in the accretion curtain is a sub-Alfvonic inflow, as opposed to the 
trans- Alfvenic outflow associated with a magnetically influenced wind. The poloidal and toroidal components of the induction 
equation yield 

Vp = /tBp, (2) 

0-^=a, (3) 
w 

where k is a function of position, Q is the angular velocity of material and a is the angular velocity of a poloidal field line. 
Since the poloidal field originates in the degenerate accretor, where it is frozen to the highly conducting stellar material, the 
constant a is also the stellar angular velocity O* . Conservation of mass and magnetic flux lead to the quantity 

e = pK=ijf, (4) 

being conserved along poloidal field-streamlines. It is noted that, since Vp and Bp are anti-parallel in the accretion curtain, 
for z > 0, K and e are negative. 

Another quantity conserved along field-streamlines follows from the angular momentum equation as 

= (5) 

This gives the total rate of transport of angular momentum, carried jointly by the gas and magnetic stresses. Equations 
(3)-(5) combine to yield 

1 - noe^/p 



and 

_ lioP/ru — fj,o\e\QtVj 
l-/xoeVp ■ 

It follows from (4) that 



(7) 



.2 „,,2 „,2 



Moe _ pVp 



(8) 



P B'^/fio vi 

where Va is the Alfven speed. The dimensionless quantity Vp/v% gives the ratio of the poloidal kinetic energy density to the 
poloidal magnetic energy density. This will be shown to be of fundamental importance to the properties of the curtain flow. 

The flow is taken to be isothermal, with sound speed a. Consideration of the rate of work done by the forces acting on 
the gas leads to the generalized Bernoulli integral 

ivl + ^^'n'-j-^^^^^ + a'lnp-n,^'n = E. (9) 



4 a G. Campbell 



This represents the sum of the poloidal and toroidal kinetic energies, the gravitational energy, and the work done by the 
pressure gradient and the magnetic torque, all per unit mass. The total energy E is conserved along poloidal magnetic field 
lines on which z = z{w). 



3 CONDITIONS AT THE CURTAIN FLOW BASE 



3.1 The poloidal flow speed 

The disc and curtain flows must be matched in the region of the disruption radius. It was shown in Campbell (2010) that the 
disc becomes disrupted over a narrow region, of width 5, inside a radius of Wra given by 

2/7 

Wm=\—\ (1-C''")'^'WA, (10) 



where 7 and em are dimensionless quantities related to the disc vertical shear distribution and difFusivity, while zua is a 
spherical Alfven radius, with 



2/7 



^ """o-"- . . (11) 

Values of 7/em ~ 1 ensure \B^/Bz\ < 1 at the base of the curtain flow, this being consistent with strong field channelling in 
a sub-Alfvenic flow. The stellar rotation parameter ^ is 



Ok 



where the corotation radius is 



with P = 27r/f2*, fiKm = flKi'Cdm) and 



The disc disruption conditions lead to 



with 



(GMP^ 



2/3 



1/3 



"i{(C^+4C)^/^-c}] 
7/2 



2/3 



C = 2 



7 



(12) 



(13) 



(14) 



(15) 



(16) 



Cm \ ^co / 

Disruption occurs inside ujco (i.e. C < 1) ^^nd is due to the rapidly growing central plane pressure caused by magnetic heating 
increasing with decreasing distance from the star. The disc expands vertically and the accretion flow is fed into the magnetically 
controlled curtain through its base, which is a ring of width 5 situated just inside Wm, with 5 <^ Wm- Hence this thin transition 
region has a radial extent of vj-m — 5 < < vo^a and a fleld line through its centre, with equation z = zivj'), is taken to have 
the base coordinates w = wo, z = z{wo) = ho- To flrst order in ho/wo, this fleld line is vertical in the disc and hence the 
coordinate wq at which it cuts the disc surface can be taken to be the same as the coordinate at which it cuts the disc 
mid-plane z = Q, for most purposes. 

As an illustration, consider 7/em = 1 and typical neutron star X-ray binary parameters M = 1.4 M©, R = 10* m, 
M = 2 X 10"^ Mq yr-\ So = 2 x 10^ T, and P = 3 s. Equations (10)-(16) then lead to ^ = 0.8, Wm = 0.8 Wco, t^A = 0.9 t^co, 
with Toco ~ 3.5 X 10® m. Fig. 1 shows a general schematic illustration of such a system. 

At disruption, the flow divides evenly between the upper and lower halves of the disc, so only the upper half need be 
considered. Matter is fed in to the upper curtain through its base at a rate 

M 

— = 2irWjn{pVz)o5, (17) 



with the poloidal field being nearly vertical here, so Vp = Vz- The ratio (Wp/wa)o can be written as 

..2 \ „ /„ \ 7/2 



(18) 



where the last expression follows from using (lb) for Bz and (17) to eliminate {pVz)o, then employing (11) for vja to eliminate 
M. This expression can be written as 

,.2 \ / \ 7/2 
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sonic point 




TUsn ^ t*7m CTco 

[P] 

Figure 1. Schematic picture of a typical system. The vertical cross-section shows the relative positions of the disc corotation radius, 

vuco, the disruption radius, vum, and the sonic point radius, zusn, in the curtain flow. Disruption occurs over a narrow annulus of width 
S through which matter flows from the disc into the magnetically channelled accretion curtain. 



It was shown in Campbell (2010) that rom/roA < 1 and ho < d ^ vj^. As the disc expands vertically in the disruption region 
hydrostatic equilibrium starts to break down and matter enters the accretion curtain with a subsonic speed Vzo < a. The 
nearly vertical force balance in the disc, together with a < Cs due to the surface temperature being less than the central plane 

temperature, gives 

^—<1. (20) 

It therefore follows from (19) that the initial curtain flow is well sub-Alfvenic. 

Evaluating the conserved quantity pVp/Bp at the disc surface and the sonic point leads to 

VzO _ \BziT^o)\ Pan 



a Sp(t^sn) po 

Employing the components (la,b) for and in 

dz Bz 



(21) 



and solving the resulting differential equation, subject to the boundary condition z{pjQ) = /lo, gives the field line equation 



Using this, together with (la,b), to evaluate the field strength B^ — (B^ + B^)^''^ along a poloidal flux tube yields 

1/2 

(24) 



Bp = 2|B2;o| — 5- 



n / \ 2/3' 



4 V tjro 



This shows that Bp increases monotonically as zu decreases from the curtain base to the accretor. 

It will be shown that the density decreases from the curtain base to the sonic point, where it has a minimum. It then 
increases along the supersonic part of the flow. Hence it follows from (8) that the maximum value of the ratio fp/vl occurs 
at the sonic point and that this is related to its value at the curtain base by 



,,2 \ /„.2 



(25) 



Using (19) and (21) then gives 



7/2 

,2 



rom\ / O70 a \ ho \Bz{xX7o)\ , , 

^aJ \ ho Vko J S Bp{wsn)' 

This is signiflcantly less than unity and hence the poloidal speed is well sub-Alfvenic everywhere in the curtain flow. Distortions 

to the stellar poloidal magnetic field are consequently small and the use of its unperturbed form in the curtain flow is therefore 
justified. Equations (7), (8) and (24) show that the field winding ratio \B^/Bp\ decreases monotonically from the disruption 
radius to the stellar surface accretion radius. It is noted that it is the product B^Bp which is related to angular momentum 
transfer to the star, and this magnetic stress term increases with decreasing tu through the curtain flow. 
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3.2 The toroidal magnetic field and angular velocity 

In the region just inside tlie disruption radius, zum, a transition occurs from tlie liiglily diffusive disc flow to the fleld controlled 
curtain flow. In the disc, the winding of the vertical magnetic field by the vertical shear is balanced by diffusion and the 
induction equation leads to the surface toroidal field 

where x = zu/zum, n > and 7/em ^ 1 corresponds to diffusion processes maintaining moderate values of field winding (see 
Campbell 2010). In the curtain flow the creation of due to the variation of f2 along Bp is balanced by its advection due 
to the poloidal flow, yielding (7) for B^. 

At TO = Too = TOm it follows that x = 1 and (27) gives 

§7) ^^{i-e/'y (28) 

The angular velocity of material must vary from a value of f2K(wm) in the disc mid-plane a = to a value consistent with 
the curtain flow solution (3) at the base z = ho, with the magnetic field being continuous acrosss the boundary. Equations 
(2) and (3) give 

VzO ( B^is 



fio = n(TOo,0o) = n. -— . (29) 

Noting from (12) that 0* = ^^/^^Km, and using (28) for (S^s/|S^|)o, (29) yields 

no = /n. (30) 

where 

/ = 1-^^MP- (31) 

em fKm C ' 

Since the last term in (31) is small, / is slightly less than unity. It is noted that this vertical transition in Q is necessary 
since otherwise B^ would attain very large negative values at the base of the curtain flow. Such values would be inconsistent 
with strong magnetic channelling in the accretion curtain and would be unsustainable, being removed by instabilities and 
reconnection. 



4 ANGULAR MOMENTUM TRANSFER 

4.1 Angular momentum transport in the curtain flow 

The total rate of transport of angular momentum, per unit poloidal magnetic flux, is given by (5). Since /3 is conserved along 
poloidal magnetic field lines, it can be expressed in terms of quantities evaluated at the accretion curtain base (toq, ho). Noting 
that S^(too,/io) = B^{izro), (5) yields 

= lelTOofio + —B^si^o). (32) 
Ho 

Since e = pVz/Bz, the ratio of the magnetic term to the material term in (32) can be expressed as 

TOog,»s(^o) ^ fvjvz B^s \ ^ J_ [ f^y^^^^ (-33^ 

fio\e\vjlno \v^V^\Bz\J^ TOm MS.lJo ' ^ ' 

where (18) has been used to eliminate {va/Vp)o- Hence the last term can be dropped in (32) to give 

P = IR I ^0"0 = T- CIO I . (34) 

where (17) has been employed to eliminate {f>Vz)o- It follows that /3 > and this will lead to a positive accretion torque on 
the star. 

It is noted that n must change from f^Km = nK(TOm) at z = to fio at the curtain base z = h{wo)- Since fio < f^Km, 
the transition from the central plane of the disc to the accretion curtain base involves a loss of angular momentum per unit 
mass of ti7^,j(nKm — f^o). The flow in this region changes from horizontal at z = to vertical at 2: = /i(too) = ho- The angular 
momentum equation is given by 

^{^pv^r^'n) + ±{wpvz^'n) = A (p^to^I^) + ^ (^-^^'B,Bz^ , (35) 

relating the divergence of the angular momentum flux to the viscous and magnetic torques. Integrating over the region 
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< z < ho, noting that and Vz vanish at 2 = 0, yields 

1 9„ ^ d f C*"" dQ 



1 d f f ' do, f ° \ 

zu(pVz)oU7^iio w^Ba,sBzs = -r- I / pi^-^—dz — vj^ I pv^Q.dz 1 

P-o dw V Jo dw Jq J 



(36) 



This transition region of the disc has a radial extent of zum — S < zu < zum- The disc completely ends at its inner edge and 
matter is field dominated here. Hence appropriate inner edge conditions are 

p{w^-6,z) = and (|^] =0. (37a,b) 

' Wm — <5 



At the outer edge 



_J = (38) 



Integrating (36) over this region, applying (37a,b) and (38), then gives 



t^m(pVz)o'Ci7mf^O "K^m-B^s-BzO 



5 = -^TOmf^Km (i^S)^^ + ^TO^nKm, (39) 



Mo 
where 

E = 2 / pdz. (40) 

and 



M = -47r / xupv^dz (41) 
Jo 

with M being the mass transfer rate into the outer surface at tu = tUm- It was shown in Campbell (2010) that the boundary 

condition 

(-S)„.='/^ (42) 

is consistent with disc disruption due to magnetic heating, for q < 0.4. Using this, together with (17) and (33), in (39) leads 
to 

Mtul^i^Km - fio) = gMro^nKm (43) 

Noting that fio = f^*, with / close to unity, and = ^^^^fiKm, (43) gives the angular momentum balance condition 

q = l-e^^. (44) 

This can be satisfied by the values of ^ of most interest, which are nearer to 1 than to 0, since in many observed systems the 
star has been spun up to a state in which ^ is well above zero. Small values of ^ may lead to a break down in steady state 
conditions. 



4.2 The accretion torque 



The accretion torque on the star can be evaluated from 

Ta= / ^Bp-dS, 

Js 



(45) 



where the area S incorporates the regions where the upper and lower curtain flows meet the stellar surface as ring accretion 
columns, and the conserved quantity j3 is given by (34). Since the accretion curtains are symmetric about the mid-plane z = 0, 
the total accretion torque will be twice that exerted through the upper curtain flow. The surface field ratio Bss/Brs, evaluated 
where the narrow accretion column meets the star, is given by 

1/21 



Bqs 

Brs 



itan( 



itan 



R 



(46) 



where 6^ is the angle a tangent to Bps makes to the vertical and 7? is the stellar radius. For white dwarfs and neutron stars, 
RIwttl <C 1 so is small and (46) yields Brs ^ B^s- Hence the accretion column is very nearly normal to the stellar surface, 
and a unit normal vector to its base is nearly parallel to z. It follows that the magnetic flux through the accretion column 

base is 

V'm = 27r-n7,5,B^,, (47) 

where J* is the width of the circular band of cylindrical radius forming the column base. So, allowing for the contribution 
from the lower accretion column band, the integral in (45) gives 

Ta = 2^'!/'m = 47rtJ7*<5«Bz*/3. (48) 
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Then, using conservation of magnetic flux through the accretion curtain to give ix7«(5*Bz» = Wm5|-Bz|o, (48) yields 

Ta = 47rti7^<5|B,lo/? = MwgOo, (49) 

where the last equality follows from using (34) for /3. Since / is very close to unity, (30) gives Oo = f2* and hence the accretion 
torque is 

Ta = MwqO.^ z, (50) 
To first order in h{wQ)/wo, zuo can be replaced by ccm so, since n, =C^''2nK(wn,), (50) can be expressed as 

Ta = C^'''Mm^nK(cc'm)z. (51) 

Noting that ^ < 1, this expression is a factor ^^^^ smaller than the standard form previously adopted for the accretion torque. 
The difference arises from the specific angular momentum 'Ci7^(f2K — H*) transferred to the disc via viscous stresses as material 
flows into the accretion curtain. For values of ^ in the lower part of the range 0.7 < ^ < 1, Ta can be significantly reduced. 



5 THE SLOW MAGNETOSONIC POINT 
5.1 The slow point speed 

The slow magnetosonic point in the curtain flow, and its relation to the sonic point, is now considered. The energy integral 
given by (9) is conserved along field lines on which z = z{w), so substitution of (4) and (6) for Vp and f2 in (9) leads to the 
relation 

H{w,p) = E. (52) 

Using dH = leads to 

dp ^ dH/dw 

dw dH/dp ■ ^ ' 

The condition dH/dp = generally yields Vp equal to the slow and fast magnetosonic speeds (e.g. Campbell 1997, Mestel 
1999). In the case considered here the flow is sub-Alfvenic, so only the slow speed solution Vp = Vs\ will result. At this point 

dH/dw must vanish, and p has a minimum value here. 
The function H(vj,p) can be expressed as 

Since the flow is well sub-Alfvenic, distortions of Bp from its unperturbed form will be small. Equation (6) gives 



Using this in (54) to form dH/dp = 0, and then (4) to eliminate Bp, leads to 



(55) 



V Psl J (1 - P-Of^/pAf 

Employing (4), (8), (30) with / = 1, and (34) yields 

vll ^ O^rogi / VpV [(rop/rosi)^ - 1]^ 

«^ UaAi [1 - («p/«a)?,]3 ■ ^ ' 

Since (wp/wa)si li this term can be dropped in the denominator of the last term in (57). It will be shown that this term is 
small, so (57) yields the first order expression 



O.CiJgl / Vp 

Vs\ = a\l- ' 



2„2 / \ 4 r 2 n 2 



2a^ V Va 



^0 _ -I 



(58) 



Hence Vsi is slightly below the sound speed a. 



5.2 The slow point coordinates 



The condition dH/dva = determines the coordinates of the slow magnetosonic point. Equations (4), (8), (23), (30), (34), 
(54) and (55) then lead to 



2 / TO aOp \ , o2 2 



2 \ 4 / 2 \ 2 



2/3 

TOO \ 2 



3 V^si 



<o = 0. (59) 
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Since (ro|dBp/dro|/Bp)ai ~ 1, (Vp/wDsi < 1, vJi/n^ro^i < 1 and ^si/wko < li (59) reduces to 



2 2 2 / roo \ 

SZ^TOsi — - I I 

6 



2/3 



Vko = 0. 



Noting that, to first order in ho/'OJo, t^o = 'cOm together with f2* = ^^''^Okiii, (60) leads to 



3/8 



and then the field line equation (23) yields 



2sl = 



1/4 



^3/4 



1/4 



TOO- 



^3/4 + ^2 



1/2 



(60) 



(61) 



(62) 



Equations (61) and (62) give the coordinates of the slow magnetosonic point. 

The condition zu^i < zuq is necessary to ensure that the slow point lies outside the disc, in the curtain fiow. Applying this 
to (61) leads to a critical value for ^ given by 



1/3 



(63) 



The slow point will lie above the disc for ^ > fcrit- The term in /iq/too niust be retained in (62) to ensure that Zsi —^has 

C Ccrit- 

The foregoing analysis has shown that for ^ ^ ^crit the slow magnetosonic point merges with the disc surface. The 
consequences of this result are discussed below. 



5.3 The slow magnetosonic and sonic points 

The quantity Vp /vj^ gives the ratio of the poloidal kinetic energy density to the poloidal magnetic energy density. This ratio is 
significantly less than unity throughout the curtain flow and hence distortions of the stellar magnetic field are small, so Bp is 
close to its unperturbed form. Since Vp/v% <C 1, (8) shows that the p-dependent term in the denominator of (7) is negligible 
and hence 

B4, = ^ ~ lj,o\e\n.zu. (64) 

This expression, together with (24) for Bp, shows that B^/Bp decreases monotonically from the curtain flow base to the star. 
Since {B^/Bp)o < 1, it follows that field winding is small as expected for a well sub-Alfvenic flow. The function 0(w, p) is 
obtained from (6), (8), (30) and (31) as 

7 v^o (l-e^/2)rooPo" 



(65) 



Em Vkui C^/^ 

Since Vzo/vKm <C 1, and zuopo/w^p decreases with inside the sonic point, the p-dependent term in (65) is always small. 
This corresponds to inertia being small relative to magnetic force and so material is effectively channelled by the magnetic 

field and nearly corotates with the star, with Q, approaching Cl, from below as vjqPo/zu'^p decreases with -cu. This corresponds 
to f2 being weakly dependent on p so, ignoring its p dependence, (54) leads to the critical point condition 

dp p^ p p p 

yielding Vp = a. This corresponds to the correction term in (58) being ignorable, due to the smallness of {Vp/v^)si, with the 
slow magnetosonic and sonic points merging to this degree of accuracy. 

With n close to f2*, the second term in (54) becomes the centrifugal potential which, together with the gravitational 
potential, dominates the poloidal kinetic energy term. The condition {dH / dzu) p = then becomes 

dzu V i9ro j ~^ \ dz ) dvj V dvj ) ^ B^ \dz ) ^' ^^'^^ 



with the effective potential 



^ = - (^2 + ,2)1/2 - 2"*^ ' (68) 



and z = z{zu) along the poloidal magnetic field. This is equivalent to the condition 

Bp • W = 0, (69) 

Equation (69) shows that at this critical point the components of the stellar gravitational force and centrifugal force along Bp 
cancel. For tu < Wsi gravity becomes increasingly dominant along Bp and the supersonic magnitude of Vp approaches free-fall 
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values. The foregoing analysis illustrates that strong magnetic field channelling of the flow leads to the slow magnetosonic 
and sonic points merging to high accuracy. 
The condition 

//(roo,Po) = -H"(Wsn,Psn) (70) 

can be used to determine the ratio Psn/po- Using (54) for H{zu,p), together with tUsn = tUsi and Zsn = Zb\ in (61) and (62), 
yields 



in 



Psn 
V PO 



1 w|o 

2 a2 



where 



/(O = r - 4 



(I) 



3/4 



t3/4 



+ 2. 



The term in v'^q/o^ makes a small contribution in (71) and can be ignored to a good approximation, giving 



Psn _ 
PO 

Using this in (21), and employing (24) for Bp, leads to 



: exp 



VzO 

a 



3/4 



^9/4 



1/4 



^3/4 



-1/2 




exp 





(71) 



(72) 



(73) 



(74) 



For ^ ^ ^crit, the sonic point merges with the curtain flow base and v^o — a. 

This illustrates the relevance of the results of the previous section. For stellar rotation rates having ^ well above ^crit, 
there is a signiflcant effective potential barrier that material must climb in order to flow through the sonic point and become 
gravitationally accelerated towards the stellar accretion poles. For values of $1* that correspond to ^ ^ Ccrit, this barrier is 
effectively removed due to the sonic point reaching the disc surface. This occurs because for these rotation rates the lower 
centrifugal force presents less of a barrier to the flow. 



6 DISCUSSION 
6.1 Summary 

The structure of an accretion curtain flow, formed as a result of the magnetic disruption of a disc, has been considered. 
Disruption occurs due to growing magnetic heating which causes the disc to expand vertically over a narrow radial region 
and matter is fed into the accretion curtain. The strong stellar magnetic field channels the flow, which subsequently passes 
through a sonic point beyond which it is gravitationally accelerated along the poloidal magnetic field to reach ncarlj' free-fall 
speeds before it accretes onto the star's surface. The angular momentum transferred to the star results in an accretion torque 
which causes the star to spin-up. This transfer occurs through the curtain flow via associated magnetic stresses. 

Matching the magnetic field across the boundary connecting the diffusive disc fiow and magnetically channelled flow 
allows the angular velocity distribution in the curtain flow to be calculated. The resulting sonic point coordinates are found, 
with the density reaching a minimum value at this critical point. The sonic point position depends on the stellar rotation 
rate, and a critical rate results below which the sonic point is shown to merge with the disc surface. Above this rotation rate, 
material leaving the disc must surmount a potential barrier in order to pass through the sonic point and be gravitationally 
accelerated towards the star along field lines. Pressure gradient forces allow material to reach the sonic point, beyond which 
they make an increasingly small contribution to the acceleration. As the critical rotation rate is approached very little thermal 
assistance is required. 

The accretion torque is calculated due to the transfer of angular momentum to the star via magnetic stresses operating 

through the curtain flow. It is shown that the angular momentum balance is maintained by some angular momentum being 
transferred back to the disc through the thin transition region over which disruption occurs. The enhanced radial gradient 
in the angular velocity occurring in this region facilitates the transfer, via viscous stresses. As a result, the accretion torque 
can be significantly lower than in the standard model in which all the angular momentum of material is transferred from the 
inner region of the disc to the star. 

A dipolar stellar magnetic field is considered here, since this will be the dominant multipole in neutron star and white 
dwarf accretors. However, the present work should have relevance to stars in which higher multipole components are significant, 
such as in T Tauri systems. Non-ax;isymmetric effects, such as azimuthal exchange of angular momentum, and time-dependent 
effects will tend to average out over a rotation period of the star. The stellar spin evolution will be a response to such an 
averaged torque. 
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6.2 Comparison with other work 

The above results can be compared with previous models of angular momentum transfer from the inner part of a magnetically 
disrupted disc to the accreting star. The present work agrees with the Ghosh & Lamb (1978, 1979a,b) picture of a thin 
inner radial boundary layer over which disruption occurs. This was shown in Campbell (2010) to be due to rapidly increasing 
magnetic dissipation in the disc, which is related to the steeply increasing value of Bp inside the corotation radius. The 
consequent heating of the disc causes it to expand vertically over a narrow region, with the vertical component of stellar 
gravity and the toroidal field magnetic pressure gradient not being able to balance increases in the vertical thermal pressure 
gradient. This confirms the idea of Campbell (1992) that disc disruption is a thermally related mechanism involving break 
down of the disc vertical equilibrium. A transition from the disc flow to the magnetically channelled curtain flow occurs 
through this inner region. 

The present work also agrees with the Ghosh & Lamb picture that the angular velocity of disc material changes from a 
Keplerian value to the stellar value across the inner boundary layer. Ghosh & Lamb assumed that all the angular momentum 
of material in this layer is transferred to the star via the field-channelled accretion flow, and this leads to the standard form 
of Mti7^nK(rom) for the magnitude of the accretion torque. However, it is shown here that the increased values of |dn/dro| in 
the bounday layer allow some angular momentum to be fed back into the disc. This corresponds to the angular momentum 
lost by matter as it flows from the central plane of the disc, where Cl = ClKi^m), to the base of the accretion curtain, where 
is close to n* Consequently, the accretion torque is reduced by a factor of where £, = rom/roco. This result is consistent 
with the continuity of the magnetic fleld across the curtain flow base, and it avoids the problem of a large value of \B^\ just 
above the disc. Such values of IB^I would be unsustainable and inconsistent with field channelling in the curtain flow. 

The work of Erkut & Alpar (2004) indicates that fl may approach fi* over a broader region than a boundary layer for 
values of ^ in the lower part of the range < ^ < 1. However, their analysis did not incorporate the thermal problem and 
the effects of magnetic heating. As outlined above, magnetic heating effects were shown by Campbell (2010) to lead to rapid 
disruption of the disc inside vijco, and are related to the steep increase in Bp. Nevertheless, this low ^ regime merits further 
investigation. As noted in Section 4.1, steady state conditions may break down at such low stellar rotation rates. 

The model developed by Shu et al (1994a, b) and Ostrikcr & Shu (1995) also has a reduced accretion torque. However, this 
is due to angular momentum removal via a magnetically influenced local wind developing just beyond the disruption region, 
rather than via viscous transport as in the present model. It is interesting to note that the work of Campbell (2011) illustrates 
that as the stellar spin rate approaches the equilibrium value, Poq, at which the accretion and magnetic disc torques cancel, 
a narrow region of disc thickening occurs beyond roco. This would aid the production of a local magnetically channelled wind 
and hence, for P close to Peq, an accretion curtain and wind flow may co-exist. Close to the equilibrium state values of the 
rotation parameter typically lie in the range 0.9 < ^ < 1 and so zu^ lies not far beneath zuco- As outlined in Section 1.2, the 
Shu et al model has disc disruption close to zuco, with a narrow X-type region generating a curtain- type inflow and a wind 
outflow. Hence the present model could share some features of the Shu et al model in the limit P Peq. 
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